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The  design,  experimentation,  construction,  and  simulation  of  Flapping  Micro  Air  Vehicles  are  relatively  new 
areas  of  research  being  studied  in  many  universities,  industries,  and  military  organizations.  For 
universities  in  particular,  flapping  MAVs  provide  an  area  of  research  that  is  easy  to  enter  into.  Apart  from 
the  many  grants  and  financing  offered,  this  area  of  research  is  affordable  in  both  start  up  and  ongoing 
materials  cost.  Most  of  the  materials  to  get  started  in  the  design  and  construction  of  flapping  MAVs  at  an 
entry  level  are  readily  available  and  can  be  purchased  at  a  local  hoppy  shop.  Because  of  this,  there  are 
many  emerging  studies  however;  very  few  have  the  more  advanced  capability  to  collect  data 
computationally  and  experimentally.  That  is  where  the  Air  Force  Research  Laboratory/ Air  Vehicles 
Directorate  can  excel.  They  plan  to  become  a  center  of  excellence  in  the  world  for  MAVs.  This  project  is 
a  step  in  that  direction. 

The  project  “Nonlinear  Structures  Optimization  for  Flexible  Flapping  Wing  Micro-Air  Vehicles”  was  funded 
with  Chief  Scientist  Innovative  Research  funds.  This  project  was  divided  into  4  parts,  Preliminary  design  of 
flapping  wing  device,  Couple  Abaqus  with  aerodynamic  model,  Comparison  of  experimental  and 
computational  results,  and  Evaluate  CSIRF  progress  for  transition  to  FY09  6.2  program.  These  four  parts 
are  discussed  below. 

Part  1 

Milestone  1:  Preliminary  design  of  flapping  wing  device,  measuring  wing  position  and  structural 
deformations,  measuring  6  DOF  forces  and  moments 

As  part  of  this  project,  a  Masters  student,  Lt  Craig  Svanberg,  focused  his  thesis  on  the  proof  of  concept  and 
development  of  a  flapping  wing  device,  various  testing  and  data  collection  techniques  for  these  small  scale 
flapping  wings  [1].  The  preliminary  design  of  a  flapping  wing  device  was  accomplished  in  the  4th  quarter  of 
FY07.  The  wing’s  position  and  structural  deformations  were  able  to  be  captured  via  high  speed  camera 
system;  however  no  scale  has  been  set  to  capture  the  exact  values.  Only  1  DOF,  the  normal  force  for  the 
flapping  mechanism,  was  able  to  be  accurately  captured.  The  other  5  DOFs  were  too  noisy  with 
interference,  however,  future  research  has  been  planned  to  accurately  capture  the  remaining  5  DOFs. 

This  study  has  increased  both  AFIT’s  and  AFRL’s  understanding  and  capabilities  for  future  research.  Lt 
Svanberg’s  Thesis,  “Biomimetic  Micro  Air  Vehicle  Testing  Development  and  Small  Scale  Flapping-Wing 
Analysis,”  can  be  seen  as  Attachment  A:  Report  Number  AFIT/GAE/ENY/08-M27. 

Part  2 

Milestone  2:  Couple  Abaqus  with  aerodynamic  model  as  a  computational  tool  modeling  flapping  wing  flight 
through  temporal  sub-iterations  or  through  cyclic  iterations 

Apart  from  physical  experiments  in  a  laboratory  setting,  computer  simulations  can  also  provide  a  wealth  of 
data  and  knowledge.  An  accurate  computer  simulation  can  not  only  give  one  the  results  of  an 
experimental  lab  set  up,  but  also  make  available  a  viable  method  for  optimization. 

Abaqus  is  a  high-fidelity  general-purpose  finite  element  capability  that  AFRL/RB  uses.  This  software  is 
capable  of  linear  or  non-linear,  implicit  or  explicit,  and  dynamic  or  static  solutions.  Abaqus  however  does 
not  incorporate  aerodynamic  loads  in  its  solution.  A  general  static  or  dynamic  solution  does  not  require 
aerodynamic  loads  in  order  to  be  accurate.  However,  for  a  computational  simulation  of  flapping  wing  flight 
to  be  precise,  aerodynamic  loads  will  need  to  be  present  to  solve  for  the  correct  forces,  moments, 
deformations,  lift,  drag  and  thrust  applied  to  the  modeled  wings. 

For  the  aero  code,  ultimately  a  high-fidelity  code  that  can  handle  low  Reynolds  numbers  would  be  idea. 
AFRL/RB  is  currently  working  on  low,  medium  and  high  fidelity  codes  that  can  handle  these  low  Reynolds 
numbers.  At  the  time  of  this  project,  a  fully  completed  aero  code  was  not  available.  So  while  Abaqus  was 
not  able  to  be  fully  coupled  with  an  true  aerodynamic  model,  proof  of  capability  for  Abaqus  to  be  coupled 
with  an  aero  model  was  verified  still  able  to  be  verified.  The  capability  for  Abaqus  to  be  coupled  with  an 
external  aerodynamics  code  was  examined  using  Abaqus  to  build  and  run  the  FEM  model,  an  estimated 
Aero  Code  to  predict  a  10%  resisting  load  to  the  model,  and  Python  Scripting  to  wrap  around  everything. 
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Building  the  Model  in  Abaqus  CAE 

The  flapping  wing  in  the  model  consisted  of  a  0.02  inch  diameter  rod,  2.5  inches  long  attached  to  an 
actuator  at  its  root.  The  actuator  was  set  to  a  periodic  flapping  of  10  Hz.  The  simulation  was  only  taken 
out  through  1  period  of  time,  t  =  0.1  seconds  with  a  peak  total  amplitude  of  114  degrees.  This  Abaqus/aero 
code  “solution”  was  developed  through  cyclic  iterations  rather  than  temporal  sub-iterations. 

Create  an  Input  file 

By  opening  Abaqus  CAE,  the  graphical  user  interface,  a  model  of  the  whole  system  or  job  can  be  created. 
This  includes  the  wing  structure,  mesh,  elements,  boundary  conditions,  actuations,  loads,  output  requests, 
and  time  steps.  Once  the  model  is  developed  in  Abaqus  CAE,  an  Input  file  can  be  written.  The  input  file 
can  be  opened  and  edited  in  Notepad;  WireRod.inp  can  be  seen  in  Attachment  B. 

Create  a  grid 

In  the  input  file,  the  grid  is  made  up  of  node,  x,  y,  and  z  positions.  For  the  Rod  part  it  looked  as  follows: 
*Part,  name=WireLine 


Node 

X 

Y 

Z 

1, 

0., 

0., 

0. 

2, 

0.5, 

0., 

0 

3, 

1., 

0., 

0. 

4, 

1.5, 

0., 

0 

5, 

2., 

0., 

0. 

6, 

2.5, 

0., 

0 

Create  a  kinematics  file 

Abaqus  needs  a  kinematics  equation  or  tabulated  data  entry  for  the  root  of  the  wing  to  actuate  as  desired. 
In  this  case,  tabulated  data  was  used  fora  1 -dimensional  10  Hz  sin  wave  for  one  period  t  =  0.0-0. 1  seconds. 
There  are  51  discrete  time  values  in  this  example. 

*Amplitude,  name=sin10tab  (time,  amplitude) 

0.,  0.,  0.002,  0.125333,  0.004,  0.24869,  0.006,  0.368125,  0.008,  0.481754,  0.01 , 0.587785,  0.012, 
0.684547,  0.014,  0.770513,  0.016,  0.844328,  0.018,  0.904827,  0.02,  0.951057,  0.022,  0.982287 
0.024,  0.998027,  0.026,  0.998027,  0.028,  0.982287,  0.03,  0.951057,  0.032,  0.904827,  0.034,  0.844328, 
0.036,  0.770513,  0.038,  0.684547,  0.04,  0.587785,  0.042,  0.481754,  0.044,  0.368125,  0.046,  0.24869 
0.048,  0.125333,  0.05,  1.22515e-16,  0.052,  -0.125333,  0.054,  -0.24869,  0.056,  -0.368125,  0.058, 
-0.481754,  0.06,  -0.587785,  0.062,  -0.684547,  0.064,  -0.770513,  0.066,  -0.844328,  0.068,  -0.904827,  0.07, 
-0.951057,  0.072,  -0.982287,  0.074,  -0.998027,  0.076,  -0.998027,  0.078,  -0.982287,  0.08,  -0.951057, 
0.082,  -0.904827,  0.084,  -0.844328,  0.086,  -0.770513,  0.088,  -0.684547,  0.09,  -0.587785,  0.092, 
-0.481754,  0.094,  -0.368125,  0.096,  -0.24869,  0.098,  -0.125333,  0.1,  -2.4503e-16 

Create  a  Forces-at-Nodes  file 

In  this  example  there  are  6  nodes.  Node  1  is  at  the  root  of  the  wing,  and  will  carry  no  forces.  Nodes  2-6 
will  carry  a  force  for  each  of  the  51  discrete  times,  t.  Tabulated  data  in  the  form  of  51  x  (time,  force)  for 
each  node  is  created.  Likewise,  moments  can  be  created  at  each  node  and  time. 

Node  2:  51  x  (time,  force) 

Node  3:  51  x  (time,  force) 

Node  4:  51  x  (time,  force) 

Node  5:  51  x  (time,  force) 

Node  6:  51  x  (time,  force) 

Start  Loop: 

Start  Python  Script 

A  python  script  can  hold  multiple  command  line  prompts. 

Run  Input  file 

Run  the  input  file  (WireRod.inp):  Abaqus  input=WireRod 
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Write  U,  V  and  A  to  file 

The  python  script  that  writes  one  or  all  of  displacements,  velocities,  accelerations,  rotational  displacements, 
rotational  velocities,  and  rotational  accelerations,  (U,  V,  A,  UR,  VR,  and  AR)  is  called  UVAatNodes.py  and 
can  be  seen  in  Attachment  C.  A,  UR,  VR,  and  AR  were  commented  out  in  this  case  using  the  #  symbol. 
Notice  that  only  U  and  V  were  written  to  the  file  UVA.txt  in  Attachment  D. 
tStep-1 

Node  =  1  U[x]  =  0.00000000,  U[y]  =  0.00000000,  u[z]  =  0.00000000 

Node  =  2  U[x]  =  0.00000000,  U[y]  =  0.00000000,  u[z]  =  0.00000000 

Node  =  3  U[x]  =  0.00000000,  U[y]  =  0.00000000,  u[z]  =  0.00000000 

Node  =  4  U[x]  =  0.00000000,  U[y]  =  0.00000000,  u[z]  =  0.00000000 

Node  =  5  U[x]  =  0.00000000,  U[y]  =  0.00000000,  u[z]  =  0.00000000 

Node  =  6  U[x]  =  0.00000000,  U[y]  =  0.00000000,  u[z]  =  0.00000000 

tStep-2 

Node  =  1  U[x]  =  0.00000000,  U[y]  =  0.00000000,  u[z]  =  0.00000000 
Node  =  2  U[x]  =  -0.00010812,  U[y]  =  0.01039741,  u[z]  =  0.00000000 

Node  =  3  U[x]  =  -0.00015341,  U[y]  =  0.01712761 ,  u[z]  =  0.00000000 

Node  =  4  U[x]  =  -0.0001 6988,  U[y]  =  0.021 1 8577,  u[z]  =  0.00000000 

Node  =  5  U[x]  =  -0.00017518,  U[y]  =  0.02348901 ,  u[z]  =  0.00000000 

Node  =  6  U[x]  =  -0.0001 7691 ,  U[y]  =  0.02480394,  u[z]  =  0.00000000 

Aero  code  reads  U,  V,  A 

Aero  code  writes  Forces-at-Nodes 

The  aero  code  reads  in  the  position,  velocity,  and  acceleration  of  each  node  at  each  time  step  finds  its 
solution  and  in  turn  outputs  a  file  of  Node,  Force  in  x  direction,  Force  in  y  direction,  and  Force  in  z  direction 
(Fx,  Fy,  and  Fz)  for  each  of  the  51  time  steps.  These  values  of  forces  and  moments  at  each  node  are  all 
Abaqus  will  need  to  update  the  next  iteration.  An  estimation  of  an  approximately  10%  resisting  force  was 
applied  in  the  opposite  direction  of  the  velocity  at  each  time  step  and  node.  DLOAD.txt  can  be  seen  in 
Attachment  E. 

tStep-1 

Node  =  1  F[x]  =  0.00000000,  F[y]  =  0.00000000,  F[z]  =  0.00000000 

Node  =  2  F[x]  =  0.00000000,  F[y]  =  0.00000000,  F[z]  =  0.00000000 

Node  =  3  F[x]  =  0.00000000,  F[y]  =  0.00000000,  F[z]  =  0.00000000 

Node  =  4  F[x]  =  0.00000000,  F[y]  =  0.00000000,  F[z]  =  0.00000000 

Node  =  5  F[x]  =  0.00000000,  F[y]  =  0.00000000,  F[z]  =  0.00000000 

Node  =  6  F[x]  =  0.00000000,  F[y]  =  0.00000000,  F[z]  =  0.00000000 

tStep-2 

Node  =  1  F[x]  =  0.00000000,  F[y]  =  0.00000000,  F[z]  =  0.00000000 
Node  =  2  F[x]  =  -0.00578963,  F[y]  =  0.08462580,  F[z]  =  0.00000004 

Node  =  3  F[x]  =  -0.00936150,  F[y]  =  0.15126581 ,  F[z]  =  0.00000089 

Node  =  4  F[x]  =  -0.01221042,  F[y]  =  0.21087037,  F[z]  =  0.00000190 

Node  =  5  F[x]  =  -0.01508493,  F[y]  =  0.27074319,  F[z]  =  0.00000251 

Node  =  6  F[x]  =  -0.00230485,  F[y]  =  0.03090382,  F[z]  =  -0.00000002 

Python  Script  updates  Forces-at-Nodes  values  in  the  Rod-Input  file. 

End  Loop  when  the  solution  converges  on  itself. 

Convergence  Verified 

Results  from  the  first  and  second  run  of  the  loop  for  Node  #6  are  shown  below.  After  the  1 0%  resisting 
force  is  applied  to  the  Abaqus  model  it  is  clearly  seen  that  the  solution  changes.  When  the  new  Abaqus 
solutions  are  fed  into  the  aero  code,  the  newer  forces  were  less  in  magnitude  than  the  initial  forces.  As  a 
result  the  3rd  runs  solution  in  this  simple  case  will  fall  in  between  the  previous  two  solutions.  This  example 
converged  very  well  after  only  5  or  6  runs. 
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From  this,  it  can  be  said  that  once  a  viable  aero  code  is  available,  coupling  Abaqus  with  the  code  using  the 
previous  method  will  be  attainable.  The  CPU  processing  time  for  a  true  aero  code  such  as  a  Vortex  Lattice 
method  solver  will  take  more  than  95%  of  the  processing  time  for  any  convergence. 

Part  3 

Milestone  3:  Comparison  of  experimental  results  and  computational  results  complete 

The  experimental  flapping  results  do  not  include  any  exact  position  or  deformations  to  be  compared  to 
computational  results.  Any  computational  results  would  include  position  and  deformations  at  any  node 
desired.  However  the  Abaqus  results  would  not  include  aero  loads  as  access  to  a  low-Reynolds  number 
medium  fidelity  aero  code  was  not  available  as  previously  thought.  Two  and  three  dimensional  vortex 
lattice  codes  are  currently  being  worked  on  in-house  at  AFRL/RB.  There  will  be  plenty  of  opportunities  in 
the  future  to  verify  and  validate  results  of  both  experimental  and  computational  results. 

Part  4 

Milestone  4:  Evaluate  CSIRF  progress  for  transition  to  FY09  6.2  program 

The  work  in  this  project  has  already  led  to  a  myriad  of  opportunities  for  improvement  and  transition  to  future 
6.2  work.  Studies  are  currently  ongoing  with  the  next  year  group  of  AFIT  students  to  obtain  all  6  DOF 
accurately,  capture  precise  position  and  deformation  values,  and  develop  more  advanced  flapping  wing 
devices.  Whether  or  not  testing  can  be  done  in  a  vacuum  chamber  is  also  being  looked  into.  If  this  is 
probable,  then  results  can  be  compared  directly  to  an  Abaqus  model  alone.  Also,  once  an  aero  code 
becomes  available,  it  should  be  relatively  straightforward  to  couple  the  two,  as  Abaqus  has  shown  that  it  is 
able  to  be  linked  external  system.  This  project  was  successfully  transitioned  to  the  6.2  project  “Hovering 
Flight  Sciences.” 

References 

[1]  C.E.  Svanberg,  “Biomimetic  Micro  air  Vehicle  Testing  Development  and  Small  Scale  Flapping-Wing 
Analysis,”  March  2008. 
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Attachment  B 


*Heading 

**  Job  name:  Linela  Model  name:  FlappingLine 
^Preprint,  echo=NO,  model=NO,  history=NO,  contact=NO 
*  * 

**  PARTS 
*  * 

*Part,  name=WireLine 
*Node 


1, 

o.. 

o.. 

0 

2, 

0.5, 

o.. 

0 

3, 

1.  , 

o.. 

0 

4, 

1.5, 

o.. 

0 

5, 

2., 

o.. 

0 

6, 

2.5, 

o.. 

0 

*Element,  type=B31 

1,  1,  2 

2,  2,  3 

3,  3,  4 

4,  4,  5 

5,  5,  6 

*Nset,  nset=_PickedSet2 ,  internal 
1,  2 

*Elset,  elset=_PickedSet2 ,  internal 

1, 

*Nset,  nset=_PickedSet5 ,  internal 
1,  2 

*Elset,  elset=_PickedSet5 ,  internal 

1, 

*Nset,  nset=_PickedSet6,  internal,  generate 
1,  6,  1 

*Elset,  elset=_PickedSet6,  internal,  generate 

1,  5,  1 

**  Region:  (CircleSection : Picked) ,  (Beam  Orientation : Picked) 
*Elset,  elset=_Il,  internal 

1, 

**  Section:  CircleSection  Profile:  circle 

*Beam  General  Section,  elset=_Il,  poisson  =  0.3,  density=0 . 00017, 
section=CIRC 
0.01 

0.,0.,-l. 

2e+07 ,  1 . 2e+07 
^Damping,  alpha=3 . 

**  Region:  (CircleSection : Picked) 

*Elset,  elset=_I2,  internal,  generate 

2,  5,  1 

**  Section:  CircleSection  Profile:  circle 

*Beam  General  Section,  elset=_l2,  poisson  =  0.3,  density=0 . 00017, 
section=CIRC 
0.01 

0.,0.,-l. 

2e+07 ,  1 . 2e+07 
^Damping,  alpha=3 . 
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*End  Part 
*  * 

*  * 

**  ASSEMBLY 
*  * 

^Assembly,  name=Assembly 
*  * 

*Instance,  name=WireLine-l ,  part=WireLine 
*End  Instance 
*  * 

*Element,  type=CONN3D2 
1,  ,  WireLine-1.1 

^Connector  Section,  elset=_PickedSet6 
Join,  Rotation 

*Nset,  nset=Wire-l-Set-l,  instance=WireLine-l 

1, 

*Elset,  elset=Wire-l-Set-l 
1, 

*Elset,  elset=_PickedSet6,  internal 

1, 

*Elset,  elset=_PickedSet8 ,  internal 

1, 

*Elset,  elset=_PickedSet9,  internal 

1, 

*Elset,  elset=_PickedSetl 1 ,  internal 

1, 

*Elset,  elset=_PickedSetl2,  internal 

1, 

*Elset,  elset=_PickedSetl3,  internal 

1, 

*Nset,  nset=_PickedSetl4 ,  internal,  instance=WireLine-l 

2, 

*Nset,  nset=_PickedSetl5,  internal,  instance=WireLine-l 

3, 

*Nset,  nset=_PickedSetl 6,  internal,  instance=WireLine-l 

4, 

*Nset,  nset=_PickedSetl7 ,  internal,  instance=WireLine-l 

5, 

*Nset,  nset=_PickedSetl 8 ,  internal,  instance=WireLine-l 

6, 

*End  Assembly 
*Amplitude,  name=nl 

0.,  0.,  0.002,  5 . 83239e-08,  0.004,  1.00638e-07,  0.006,  5.60611e-08 

0.008,  5 . 28955e-08,  0.01,  7.8951e-08,  0.012,  3.42409e-08,  0.014,  6.64435e- 

09 

0.016,  3 . 74796e-08,  0.018,  2.91513e-08,  0.02,  2.18218e-08,  0.022, 

2 . 59273e-08 

0.024,  2 . 47149e-08,  0.026,  -4 . 12702e-08,  0.028,  -2 . 72 973e-08 ,  0.03, 

5 . 08839e-10 

0.032,  -2 . 338e-08 ,  0.034,  -2 . 87364e-08,  0.036,  -1 .  84142e-08,  0.038,  - 
5 . 95595e-08 

0.04,  -8 . 32315e-08,  0.042,  -4 . 55654e-08,  0.044,  -4 . 71681e-08,  0.046,  - 
7.387  92e-08 
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0.048,  -4 . 5457e-0  8 ,  0.05,  -3 . 46535e-08,  0.052,  -8 . 61659e-08,  0.054, 
8.78447e-08 

0.056,  -4 . 52237e-08,  0.058,  -5 . 47121e-08,  0.06,  -6 . 18404e-08,  0.062, 
2 . 02594e-08 

0.064,  -2 . 01634e-08,  0.066,  -5 . 44 641e-08 ,  0.068,  -3 . 73712e-08,  0.07, 
1 . 7278e-08 

0.072,  -2 . 74173e-08,  0.074,  8.40318e-10,  0.076,  2.93238e-08,  0.078, 
6. 1483e-10 

0.08,  -6. 60815e-09,  0.082,  3.32018e-08,  0.084,  3.76589e-08,  0.086, 

3 . 251 97e-08 

0.088,  6 . 67  8  8  9e-0  8 ,  0.09,  6.71171e-08,  0.092,  2.93088e-08,  0.094, 

4 . 49846e-08 

0.096,  7.1587  6e-0  8 ,  0.098,  5.1235e-08,  0.1,  5.83141e-08 
*Amplitude,  name=n2 

0.,  0.,  0.002,  1 . 2  65e-07 ,  0.004,  1.98816e-07,  0.006,  7.85534e-08 
0.008,  1 . 18093e-07,  0.01,  2.02067e-07,  0.012,  5.16148e-08,  0.014,  - 

2 . 48024e-08 

0.016,  9 . 56447e-08,  0.018,  7.13413e-08,  0.02,  -1.3886e-09,  0.022, 

9. 16473e-08 

0.024,  7 . 56824e-08,  0.026,  -1 . 16451e-07,  0.028,  -5 . 65906e-08,  0.03, 

3 . 93204e-08 

0.032,  -6 . 17116e-08,  0.034,  -7 . 75568e-08,  0.036,  -2 . 59821e-09,  0.038 
1 . 19512e-07 

0.04,  -2 . 0557e-07 ,  0.042,  -7.7836e-08,  0.044,  -7 . 32229e-08,  0.046,  - 

1 .71753e-07 

0.048,  -8.4  637  6e-08 ,  0.05,  -4 . 102 94e-08,  0.052,  -1 . 95801e-07,  0.054, 

1 . 93055e-07 

0.056,  -6 . 53403e-08,  0.058,  -1 . 19391e-07,  0.06,  -1 . 45296e-07,  0.062, 

7 . 62328e-09 

0.064,  -2 . 82e-08 ,  0.066,  -1 . 44708e-07 ,  0.068,  -6 . 79898e-08,  0.07,  - 

1 . 10964e-08 

0.072,  -8 . 35958e-08,  0.074,  -6 . 40802e-09,  0.076,  9.25917e-08,  0.078, 

1 . 52145e-08 

0.08,  -4 . 45598e-08,  0.082,  9.12338e-08,  0.084,  8.48347e-08,  0.086, 

3 . 71343e-08 

0.088,  1 . 50064e-07,  0.09,  1.60074e-07,  0.092,  2.94592e-08,  0.094, 

7 . 91209e-08 

0.096,  1 . 68906e-07,  0.098,  8.85805e-08,  0.1,  1.01554e-07 
*Amplitude,  name=n3 

0.,  0.,  0.002,  2 . 21858e-07,  0.004,  3.09074e-07,  0.006,  6.46338e-08 
0.008,  1 . 6045e-07 ,  0.01,  3.24999e-07,  0.012,  5.69442e-08,  0.014,  - 
2 . 39885e-08 

0.016,  2 . 10052e-07,  0.018,  1.00388e-07,  0.02,  -8 . 01637e-08,  0.022, 

1 .55226e-07 

0.024,  1 . 31244e-07,  0.026,  -1 . 82612e-07,  0.028,  -6 . 06676e-08,  0.03, 

1 . 01833e-07 

0.032,  -1 . 35001e-07 ,  0.034,  -1 . 68979e-07,  0.036,  3.12797e-08,  0.038, 

1 .55522e-07 

0.04,  -3.27  67  6e-0  7 ,  0.042,  -9 . 31111e-08,  0.044,  -8 . 78 651e-08 ,  0.046, 
3 . 004  99e-07 

0.048,  -1 ,45466e-07,  0.05,  -4 . 37745e-08,  0.052,  -3 . 04522e-07,  0.054, 

2 . 79682e-07 

0.056,  -5 . 43113e-08,  0.058,  -1 . 94451e-07,  0.06,  -2 . 55138e-07,  0.062, 

5 . 38153e-09 
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0.064,  -4 . 92144e-08,  0.066,  -2 . 53 992e-07 ,  0.068,  -6 . 52349e-08,  0.07, 

3 . 49325e-08 

0.072,  -1 ,57326e-07,  0.074,  -3 . 56408e-08,  0.076,  1.67005e-07,  0.078,  - 
3.27  62  6e-08 

0.08,  -8 . 18388e-08,  0.082,  1.82141e-07,  0.084,  1.36893e-07,  0.086, 

2 . 04098e-09 

0.088,  2 . 2  008  6e-0  7 ,  0.09,  2.63886e-07,  0.092,  2.55907e-08,  0.094,  1.2675e 

07 

0.096,  2 . 98441e-07,  0.098,  1.18898e-07,  0.1,  1.19129e-07 

*Amplitude,  name=n4 

0.,  0.,  0.002,  3 . 25901e-07,  0.004,  4.36289e-07,  0.006,  3.29171e-08 
0.008,  1 . 74755e-07,  0.01,  4.42952e-07,  0.012,  6.34649e-08,  0.014,  - 

1 .59914e-08 

0.016,  3 . 66981e-07,  0.018,  1.31274e-07,  0.02,  -1 . 697 68e-07 ,  0.022, 

1 . 72577e-07 

0.024,  1 . 75828e-07,  0.026,  -2.3649e-07,  0.028,  -5 .21735e-08,  0.03, 

1 . 61458e-07 


0.032,  -2 ,39518e-07, 
1 . 7  98  8e-07 

0 

.034, 

-2 . 85247e-07, 

0.036 

,  7 . 04064e-08, 

0.038,  - 

0.04,  -4 . 31022e-07, 
4.71 8  05e- 07 

0  . 

042, 

-7 . 66451e-08, 

0.044, 

-1 . 03855e-07, 

0.046,  - 

0.048,  -2 ,28747e-07, 
3 . 50078e-07 

0 

.05, 

-3 . 88744e-08, 

0.052, 

-4 . 11 647e-07 , 

0.054,  - 

0.056,  -4 . 63906e-09, 
4 . 9351e-09 

0 

.058, 

-2 . 65  654e-07 , 

0.06, 

-4.006  9  6e- 07 , 

0.062,  - 

0.064,  -7 .582  66e- 08 , 
1 . 11461e-07 

0 

.066, 

-3 . 64  68  9e-07 , 

0.068 

,  -4 . 02302e-08 

l> 

o 

o 

0.072,  -2 . 31552e-07, 
3 . 72107e-08 

0 

.074, 

-7 ,74686e-08. 

0.076 

,  2 . 40842e-07, 

0.078,  - 

0.08,  -9. 98661e-08, 

6 . 40358e-08 

0  . 

082, 

2 . 96739e-07,  0 

.  084, 

1 . 84  68  6e-07 ,  0 

.086,  - 

0.088,  2 . 73873e-07, 
07 

0  . 

09,  3 

. 64889e-07,  0. 

092,  2 

. 30882e-08 ,  0. 

094,  2 . 0002e 

0.096,  4 . 60319e-07, 

0  . 

098, 

1 . 41243e-07,  0 

\ — 1 

\ — 1 

07836e-07 

*Amplitude,  name=n5 

0.,  0.,  0.002,  4 . 71521e-07,  0.004,  5.89313e-07,  0.006,  -5.84963e-08 
0.008,  1 . 23528e-07,  0.01,  5.91595e-07,  0.012,  1.28318e-07,  0.014,  - 
1 .34961e-08 

0.016,  4 . 95595e-07,  0.018,  2.00777e-07,  0.02,  -2 . 65703e-07,  0.022, 


1 . 75302e-07 

0.024,  2 . 07432e-07,  0.026, 

-2 . 87356e-07, 

0.028, 

-4 . 96732e-08, 

0.03, 

2 . 07247e-07 

0.032,  -3 . 59772e-07,  0.034, 

-4 . 03405e-07, 

0.036, 

1 . 10279e-07, 

0.038 

2 . 08935e-07 

0.04,  -5.20654e-07,  0.042, 

-3 . 16689e-08, 

0.044, 

-1 . 24  04 le-07 , 

0.046 

6 . 7539e-07 

0.048,  -3 .23593e-07,  0.05, 

-2 . 32369e-08, 

0.052, 

-5.20613e-07, 

0.054 

4 . 1615e-07 

0.056,  7 . 284e-08,  0.058,  -3 .28414e-07,  0.06,  -5 . 74307e-07 ,  0.062,  - 
3 . 59483e-08 

0.064,  -9 . 62563e-08,  0.066,  -4 . 70248e-07,  0.068,  -1 .25438e-08,  0.07, 

2 . 00936e-07 

0.072,  -2 . 98303e-07,  0.074,  -1.1934e-07,  0.076,  3.12152e-07,  0.078,  - 
2 . 83845e-08 
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0.08,  -1 . 01626e-07,  0.082,  4.20262e-07,  0.084,  2.25662e-07, 
1 .43304e-07 

0.088,  3 . 21114e-07,  0.09,  4.59739e-07,  0.092,  1.98534e-08, 

2 . 917  69e-07 

0.096,  6 . 4188e-07 ,  0.098,  1.56801e-07,  0.1,  7.82015e-08 
*Amplitude,  name=sinlO,  def inition=PERIODIC 
1,  62.8319,  0.,  0. 

0.,  1. 

* Amplitude,  name=sinl Otab 


o 

o 

t 

0.002,  0 

.125333,  0.004,  0. 

24869, 

0.006,  0. 

368125 

0.008, 

0 

.481754, 

0.01,  0.587785,  0 

.012, 

0 . 684547, 

0.014, 

0. 

0.016, 

0 

. 844328, 

0.018,  0.904827, 

0.02, 

0 . 951057, 

0.022, 

0. 

0.024, 

0 

. 998027, 

0.026,  0.998027, 

0.028, 

0.982287, 

0.03, 

0. 

0.032, 

0 

. 904827, 

0.034,  0.844328, 

0.036, 

0.770513, 

0.038 

,  o 

0.04, 

0. 

587785, 

0.042,  0.481754,  0 

.044, 

0.368125, 

0.046, 

0. 

0.048, 

0 

. 125333, 

0.05,  1 . 22515e-l 6 

,  0.052,  -0.125333,  0. 

054 

0.056,  -0.368125,  0.058,  -0.481754,  0.06,  -0.587785,  0.062, 
0.064,  -0.770513,  0.066,  -0.844328,  0.068,  -0.904827,  0.07, 
0.072,  -0.982287,  0.074,  -0.998027,  0.076,  -0.998027,  0.078 
0.08,  -0.951057,  0.082,  -0.904827,  0.084,  -0.844328,  0.086, 
0.088,  -0.684547,  0.09,  -0.587785,  0.092,  -0.481754,  0.094, 
0.096,  -0.24869,  0.098,  -0.125333,  0.1,  -2.4503e-16 


**  STEP:  Step-1 
*  * 

*Step,  name=Step-l,  nlgeom=YES,  inc=100000 
* Dynamic, alpha=-0 . 05,  haftol=0 . 1 
2e-06, 0 . 1, le-07, 0 . 002 
*  * 

**  BOUNDARY  CONDITIONS 
*  * 

**  Name:  BC-1  Type:  Connector  displacement 
^Connector  Motion 
_PickedSetl 1 ,  4 

**  Name:  BC-2  Type:  Connector  displacement 
*Connector  Motion 
_PickedSetl2 ,  5 

**  Name:  BC-3  Type:  Connector  displacement 
^Connector  Motion,  amplitude=sinl Otab 

_PickedSetl3,  6,  0.2 
*  * 

**  LOADS 
*  * 

**  Name:  LI  Type:  Concentrated  force 
*Cload,  amplitude=nl 
_PickedSetl4 ,  2,  -150 

**  Name:  L2  Type:  Concentrated  force 
*Cload,  amplitude=n2 
_PickedSetl5,  2,  -150 

**  Name:  L3  Type:  Concentrated  force 
*Cload,  amplitude=n3 
_PickedSetl 6,  2,  -150 

**  Name:  L4  Type:  Concentrated  force 


0.086,  - 
0.094, 


770513 
982287 
951057 
. 684547 
24869 

,  -0.24869 
-0 . 684547 
-0 . 951057 
,  -0.982287 
-0.770513 
-0.368125 
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*Cload,  amplitude=n4 
_PickedSetl7 ,  2,  -150 

**  Name:  L5  Type:  Concentrated  force 
*Cload,  amplitude=n5 
_PickedSetl8,  2,  -150 
*  * 

**  OUTPUT  REQUESTS 
*  * 

*Restart,  write,  frequency=0 
*  * 

**  FIELD  OUTPUT:  F-Output-1 
*  * 

*Output,  field,  time  interval=0 . 002,  time  marks=NO 
*Node  Output 
A,  U,  V 
*  * 

**  HISTORY  OUTPUT:  H-Output-1 
*  * 

*Output,  history,  variable=PRESELECT ,  time  interval=0 . 005 ,  time  marks=NO 
*End  Step 
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Attachment  C 


from  odbAccess  import  * 
from  abaqusConstants  import  * 

odb  =  openOdb (path= ' 1 . odb ' ) 

filename  =  "UVA.txt" 

print  "writing  to  file  %s"  %  filename 
file  =  open ( filename,  'w') 

lastFrame  =  odb . steps [' Step-1 ']• frames [-1 ] 

for  fieldName  in  lastFrame . fieldOutputs . keys () : 
file . write ( fieldName  +  "\n") 

totaltime  =  . 1 
totalsteps  =  50 

itertime  =  totaltime/totalsteps 

################## 

#U  displacement  at  nodes  1-6 

for  i  in  range (51) : 

lastFrame  =  odb. steps [' Step-1 ']• frames [i] 

iteration  =  "iteration  %3.2d"  %(i) 
for  stepName  in  odb . steps . keys () : 

#  file .write (stepName  +  "\n") 

#  file .write (iteration  +  "\n") 

#  timex  =  itertime*i 

#  time  =  "Time  =  %6.6f"  %  (timex) 
f ile . write ( " \n" ) 

#  file . write (time  +  "\n") 

displacement=lastFrame . fieldOutputs [ ' U ' ] 
f ieldValues=displacement . values 

#  For  each  U  value,  print  the  nodeLabel 

#  and  data  members . 

for  v  in  fieldValues: 

outline  =  "Node  =  %2 . Od  U[x]  =  %6.8f,  U[y]  =  %6.8f,  u[z]  =  %6.8f 
"  %  (v . nodeLabel , 

v.data[0],  v.data[l],  v. data [2]) 
f ile . write (outline) 


################## 

#UR  displacement  rotation  at  nodes  1-6 
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# 

#  displacementrotation=lastFrame . fieldOutputs [ ' UR' ] 

#  f ieldValuesur=displacementrotation .values 

# 

#  For  each  UR  value,  print  the  nodeLabel 

#  and  data  members . 

#  for  v  in  f ieldValuesur : 

#  outline  =  "Node  =  %2 . Od  UR[x]  =  %6.8f,  UR[y]  =  %6.8f,  uR[z]  = 

%6.8f  "  %  (v. nodeLabel, 

#  v.data[0],  v.data[l],  v. data [2]) 

#  file . write (outline) 

################## 

#V  velocity  at  nodes  1-6 

for  i  in  range (51) : 

#  timex  =  itertime*i 

#  time  =  "Time  =  %6.6f"  %  (timex) 

f ile . write ( " \n" ) 

#  file . write (time  +  "\n") 

lastFrame  =  odb . steps [' Step-1 ']. frames [i ] 
velocity=last Frame . fieldOutputs [ ' V ' ] 
f ieldValuesv=velocity .values 

#  For  each  V  value,  print  the  nodeLabel 

#  and  data  members . 

for  v  in  f ieldValuesv : 

outline  =  "Node  =  %2 . Od  V[x]  =  %6.8f,  V[y]  =  %6.8f,  V[z]  =  %6.8f 
"  %  (v . nodeLabel , 

v.data[0],  v.data[l],  v. data [2]) 
f ile . write (outline) 


################## 

#VR  velocity  rotation  at  nodes  1-6 

#  velocityrotation=last Frame . fieldOutputs [ ' VR ' ] 

#  f ieldValuesvr=velocityrotation .values 

#  For  each  VR  value,  print  the  nodeLabel 

#  and  data  members . 

#  for  v  in  f ieldValuesvr : 

#  outline  =  "Node  =  %2 . Od  VR[x]  =  %6.8f,  VR[y]  =  %6.8f,  VR[z]  = 
%6.8f"  %  (v . nodeLabel , 

#  v.data[0],  v.data[l],  v. data [2]) 

#  file . write (outline  +"\n") 
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################## 

#A  acceleration  at  nodes  1-6 


#  acceleration=last Frame . f ieldOutputs [ ' A ' ] 

#  f ieldValuesa=acceleration .values 


#  For  each  A  value,  print  the  nodeLabel 

#  and  data  members . 

#  for  v  in  f ieldValuesa : 

#  outline  =  "Node  =  %2 . Od  A[x]  =  %6.8f,  A[y]  =  %6.8f,  A[z]  = 
%6.8f"  %  (v . nodeLabel , 

#  v.data[0],  v.data[l],  v. data [2]) 

#  file . write (outline  +"\n") 


################## 

#AR  acceleration  rotation  at  nodes  1-6 

#  accelerationrotation=last Frame . f ieldOutputs [ ' AR ' ] 

#  f ieldValuesar=accelerationrotation .values 


#  For  each  ar  value,  print  the  nodeLabel 

#  and  data  members . 

#  for  v  in  f ieldValuesar : 

#  outline  =  "Node  =  %2 . Od  AR[x]  =  %6.8f,  AR[y]  =  %6.8f,  AR[z] 
%6.8f"  %  (v . nodeLabel , 

#  v.data[0],  v.data[l],  v. data [2]) 

#  file . write (outline  +"\n") 


file. close () 

print  "done!  open  results  in  excel" 
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